Discontinuous liquid rise in capillaries with nonuniform cross-sections 
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We consider theoretically liquid rise against gravity in capillaries with height-dependent cross- 
section. For a conical capillary made from a hydrophobic surface and dipped in a liquid reservoir, 
the equilibrium liquid height depends on the cone opening angle a, the Young-Dupre contact angle 
9, the cone radius at the reservoir's level Ro and the capillary length k~ . As a is increased from 
zero, the meniscus' position changes continuously until, when a attains a critical value, the meniscus 
jumps to the bottom of the capillary. For hydrophilic surfaces the meniscus jumps to the top. The 
same liquid height discontinuuity can be achieved with electrowetting with no mechanical motion. 
Essentially the same behavior is found for two tilted surfaces. We further consider capillaries with 
periodic radius modulations, and find that there are few competing minima for the meniscus location. 
A transition from one to another can be performed by the use of electrowetting. The phenomenon 
discussed here may find uses in microfluidic applications requiring the transport small amounts of 
water "quanta" (volume< 1 nL) in a regular fashion. 

PACS numbers: 68.08.Bc, 68.15.+e,68.18.Jk 



The behavior of liquids confined by solid surfaces is 
important in areas such as microfluidics [ij, |2(, wetting 
of porous media yK the creation of hydrophobic surfaces 
|j,|5|, oil recovery [6} and water transport in plants |j|. As 
the system size is reduced, the interfacial tensions become 
increasingly important in comparison to bulk energies, 
and are essential in understanding the equilibrium states 
as well as the system dynamics. 

Wetting has been studied for liquids in contact with 
curved surfaces 0, M, 0, UJ| , wedges \V\L Il2| and cones 
lid Hjj lllj , and to pog raphicall y llq or chemically mod- 
ulated substrates [17J [la [lfj, |2fJ. However, surprises 
appear even for very simple geometries of the bounding 
surfaces. Here we focus on the rise of a liquid in capil- 
laries with nonuniform cross-sections. When a solid cap- 
illary is immersed in a bath of liquid, the height of the 
contact line above the bath level h is given by 




(a) 



h — CK 



cos 9/R 



(1) 



h 


/e\ \ 


a 


( Ro \ 



(b) 



where n^ 1 = (cr/gp) 1 / 2 is the capillary length, a is the 
liquid-gas interfacial tension, p is the liquid mass density 
(gas density neglected), g is the gravitational accelera- 
tion, and 6 is the Young-Dupre contact angle given by 
cos0 = (~fg S — Jis)/&, where j gs and ~fi s are the gas-solid 
and liquid-solid interfacial tensions, respectively [l|. The 
constant c is c = 2 for a cylindrical capillary, in which 
case R is the radius, or c = 1 for two parallel and flat 
surfaces separated at distance 2R. The liquid is sucked 
upwards if the capillary's surface is hydrophilic (9 < in), 
and is depressed downwards in the case of a hydrophobic 
surface (9 > |tt) HJ. 

Suppose now that the capillary walls are not vertical, 
but rather have some opening angle a as is illustrated in 
Fig. 1. What is the liquid rise then? One can naively 
expect that if a is small, h changes from Eq. (J) by a 



FIG. I: Schematic illustration of cone capillary or two tilted 
planes, and definitions of parameters. Two of the possible 
cases: (a) hydrophobic surface, cos 9 < 0, positive opening 
angle a and negative h. (b) hydrophilic surface, cos 9 > 0, 
negative a and positive h. 



small amount proportional to a; it is not even a-priori 
clear whether h increases or decreases. We restrict our- 
selves to narrow capillaries, where kR <C 1 is satisfied. 
In this case, as will be verified below, the height is larger 
than the radius, h 3> R, and the height variations of 
the meniscus surface are negligible compared to the total 
height. 

In mechanical equilibrium, at the contact line the 



Laplace pressure is balanced by the hydrostatic pressure 

CO 



Po 
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P - pgh 
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where Pq is the ambient pressure and r is the inverse 
curvature, and is given by r(h) = — R(h) / cos(8 + a) . We 
denote Ro as the radius at the bath level (see Fig. 1), 
and hence R(h) = Rq + htana. 



(a) 



i(h) 





\ i i i i 




0.4 


\ 

\ 
\ 


A 


0.2 


\ 

\ 
\ 


»_ 







' 


0.2 




- 


04 


_ 


_ 




II 



-4 -3 -2 



h 



i(h) 




FIG. 2: (a) Solid curved line is a plot of f(h) from Eq. £IJ an d 
solid horizontal line is cos(# + a). Their intersection occurs 
at two points hi (the meniscus location) and I12 < hi marked 
with circles. R = 0.5, c = 2 and a = 0.087 (5°). Dashed 
lines are the same, but a is twice as large, a = 0.174 (10°), 
above the critical angle. In this case there is no solution to 
Eq. 0, and the meniscus jumps down to the bottom of 
the capillary, (b) Plot of f(h) in the case of periodically- 
modulated capillary [Eq. ©]. The horizontal dashed line is 
cos 8, and the multiple intersections give the possible meniscus 
locations, c = 2, _Ro = 0.07, R m = Ro/2 and /1 = 3. 



We therefore find that the liquid rise is given by 

cos(9 + a) = f(h) (3) 

f(h) = -h(R Q + hta,na) (4) 

c 

where the dimensionless variable h = nh and Rq = kRq 
have been used. These equations reduce to the familiar 
form Eq. QJ in the limit a — > 0. Let us concentrate 



first on the case where a is positive and the surface is 
hydrophobic, cos 9 < [see Fig. 1 (a)]; the results for 
a < follow immediately. The left hand-side of Eq. © 
is then negative for small enough values of a, and the 
quadratic form of f(h) means that the two solutions h± 
and hi are negative, see Fig. 2 (a). The stable solution 
is h\ while hi < hi is unstable. 

If the opening angle a is too large, however, the 
minimum of /(/i), attained at h* = — i?o/(2tana), is 
f(h*) = — -Ro/(4ctana) > cos(8 + a), and there is no 
solution. Hence, for a given value of contact angle 8, 
the critical value of the opening angle a c is given by the 
condition f(h*) = cos(9 + a c ). As a is increased past 
a c , the meniscus "jumps" all the way to the bottom of 
the capillary; in the case of a nearly closed capillary this 
occurs at h = 2h*. 

When the surface is hydrophilic and both 9 and a are 
small, then there is always a positive solution for h. How- 



a > ^77, the liquid height is 



ever, if 6 < \tt but 9 

negative and the jump again is possible. In essence the 

capillary behaves as a hydrophobic surface. 

A different approach, potentially useful in applications, 
is that of electrowetting. In the experimental setup the 
opening angle a is fixed, but the contact angle may be 
changed with an external potential V imposed on the 
walls: e = 9(V). The change to cos6> is eV 2 /{2a\ D ) ~ 
0.3V" 2 (where V is in Volt) Hi HI [gfj, and thus can be 
quite large (we took the dielectric constant of water and 
the Debye screening length Ad = 10 nm). At a fixed 
value of a, an increase in 9 lowers the liquid height until 
9 reaches 9 C given by 
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(5) 



At all 9 > 9 C the meniscus jumps again to the bottom of 
the capillary. However, if a < a* , where a* is given by 

Rl/Ac , 



sin a 



(6) 



the liquid height is a continuous function of 9 at all 9. 
The threshold angle a* is quite small; if Ro = 0.1 we find 
a* = 1.25 -10- 3 (0.07°). 

Figure 3 (a) is a phase-diagram in the a-9 plane. In the 
region marked "continuous" and for positive a, h(a, 9) 
changes continuously. Across the critical line 9 c (a) [Eq. 
|J5J|], h changes discontinuously (meniscus is at the bot- 
tom of the capillary) . Figure 3 (b) shows the height h as 
a function of 9 at fixed value of a. The meniscus height 
h decreases below zero until, at the critical value of 9, its 
height jumps from h* to the capillary bottom (at 2h* if 
the capillary is nearly closed) . Further increase of 9 does 
not change the meniscus' location. 

The liquid behavior in capillaries with negative a [Fig. 
1 (b)] follows from the symmetry of the problem: the 
transformation a — > —a and 9 — > tt — 9 leaves Eqs. © 
and 10} unchanged if h — > —h. For negative values of a, a 
decrease of 9 from large values to small ones past 9 C leads 
to a jump of the meniscus to the top of the capillary. 
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FIG. 3: (a) Phase-diagram in the opening angle-contact an- 
gle (a-9) plane. For positive a, solid line is 9 c (a) from Eq. 
©, separating two regions: below it, the meniscus height h 
changes continuously as a function of a and 9. h is positive 

if h- 

9 + a — i-7r (dashed diagonal line). At a given opening angle 
a, increase of 9 past 9 C leads to a jump down of the menis- 
cus from h — h* < to the bottom of the capillary. Above 
the critical line the meniscus is at the bottom. For all angles 
a < a* (see text) the behavior is continuous. The phase- 
diagram is symmetric with respect to a — > —a, 9 — > n — 9 
and h — > —h. Parameters are Ro — 1 and c — 2. (b) Menis- 
cus location iasa function of the contact angle 9 at fixed 
a = 0.01, c = 2 and Rq = 0.1. 



The above insight can be used to exploring differ- 
ent capillaries, and we briefly mention a capillary with 
periodic width modulations [271 l28j . namely R(h) = 
i?o + R m sux(qh), where R m is the modulation amplitude 
and q its wavenumber. We restrict ourselves to the long 
wavelength regime, where qR m < 1. In this case it can 
be shown that the governing equations replacing Eqs. © 
and 10} are 

cosfl = f(h) (7) 

f(h) = -h(R + R m sm{^h)) (8) 



FIG. 4: Schematic illustration of a capillary with sinusoidal 
modulations of the radius. 



where R m = nR m and fi = q/n. It is clear from Fig. 2 
(b) that there are multiple solutions, half of which are 
maxima and the other half minima. 

For a system prepared in a given minimum, increasing 
9 by the use of electrowetting decreases cos 9. Thus, the 
liquid location changes - h decrease. When the liquid 
height overlaps with a minimum of f(h), further increase 
of 9 leads to a jump in the liquid height to the next 
"branch" of f(h). In this way one "quantum" of liquid 
is depleted from the capillary; if 9 is decreased, at each 
step one liquid unit is sucked into the capillary. The unit 
volume can be estimated to be v ~ Rq/ 'q; for a capillary 



10 d 



we find v — 10 nL, whereas reducing the sizes to Rq = 10 
/xm and q = 10 m _1 gives v = 10~ 2 nL. 

In summary, we have shown that the liquid height in 
capillaries with nonuniform cross-sections is a discontin- 
uous function of the geometrical variables. This peculiar 
phase-transition is important for the understanding of 
liquids confined to small environments, as capillaries in 
practice rarely have uniform cross-sections. Indeed, since 
a* can be extremely small, if the surfaces of the capillary 
are super-hydrophobic, very small deviations of a around 
a = will yield discontinuous liquid heights. It may be 
beneficial to exploit the dependence of the water level 
on the contact angle (e.g.) in microfluidic applications 
where it is desired to accurately control small volumes 
of liquid. The setup of Fig. 1 could possibly be used 
as a "switch" to prevent or allow liquid flow or electrical 
current in the direction perpendicular to the plane of the 
paper, while the setup of Fig. 4 permits to "suck" known 
quantities of fluids. 

In this paper we assumed ideal surfaces and thus have 
not dealt with rough surfaces and hysteresis effects [29j. 
It would also be interesting to consider a pressurized cap- 



illary (micro-pipette), where the finite pressure difference 
between the capillary and the ambient atmosphere can 
bring about further surprises |30| . 
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